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Summary 

1. The huge changes in population sizes of Arctic-nesting geese offer a great opportunity to 
study population limitation in migratory animals. In geese, population limitation seems to have 
shifted from wintering to summering grounds. There, in the Arctic, climate is rapidly changing, 
and this may impact reproductive performance, and perhaps population size of geese, both 
directly (e.g. by changes in snow melt) or indirectly (e.g. by changes in trophic interactions). 

2. Dark-bellied brent geese (Branta bemicla bernicla L.) increased 20-fold since the 1950s. Its 
reproduction fluctuates strongly in concert with the 3-year lemming cycle. An earlier analysis, 
covering the growth period until 1988, did not find evidence for density dependence, but there- 
after the population levelled off and even decreased. The question is whether this is caused by 
changes in lemming cycles, population density or other factors like carry-over effects. 

3. Breeding success was derived from proportions of juveniles. We used an information-theo- 
retical approach to investigate which environmental factors best explained the variation in 
breeding success over nearly 50 years (1960-2008). We subsequently combined GLM predic- 
tions of breeding success with published survival estimates to project the population trajectory 
since 1991 (year of maximum population size). In this way, we separated the effects of lem- 
ming abundance and population density on population development. 

4. Breeding success was mainly dependent on lemming abundance, the onset of spring at the 
breeding grounds, and the population size of brent goose. No evidence was found for carry- 
over effects (i.e. effects of conditions at main spring staging site). Negative density depen- 
dence was operating at a population size above c. 200 000 individuals, but the levelling off of 
the population could be explained by faltering lemming cycles alone. 

5. Lemmings have long been known to affect population productivity of Arctic-nesting 
migratory birds and, more recently, possibly population dynamics of resident bird species, 
but this is the first evidence for effects of lemming abundance on population size of a migra- 
tory bird species. Why lemming cycles are faltering in the last two decades is unclear, but this 
may be associated with changes in winter climate at Taimyr Peninsula (Siberia). 

Key-words: bird migration, climate change, dark-bellied brent goose, density dependence, 
reproductive success 



Introduction 

Unlike resident species, vital rates and hence the equilib- 
rium population sizes of migratory species may be 
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affected by conditions encountered on the breeding, win- 
tering, and stopover sites (Newton 2004). For resident 
and migratory species alike, a change in population size 
may in itself affect survival and reproduction through 
negative density dependence (Newton 1998). However, as 
population limitation can occur in more than one part of 
the world, achieving an understanding of density depen- 
dence is difficult in migratory species (Newton 2004). 
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Moreover, conditions experienced in one part of the 
world can potentially affect the performance of the 
migrants in other parts (carry-over effects; Marra, 
Hobson & Holmes 1998; Norris & Taylor 2006; Harrison 
et al. 2011). A multi-factorial approach is therefore 
needed to unravel the causes of variation in vital rates of 
migratory species. 

Swan and goose species (Anserini), many of which 
breed in northern regions and winter in the temperate 
zone, seem to be well suited for such a multi-factorial 
approach to the study of vital rates and equilibrium popu- 
lation sizes of migrants. Because young of the year spend 
the winter with their parents (Cramp & Simmons 1977), 
both reproductive and annual mortality rate can be quan- 
tified on the wintering grounds. Most swan and goose 
species have increased over the past 50 years, which has 
been attributed to a higher winter survival because of 
shooting restrictions (Ebbinge 1991) as well as improved 
feeding conditions related to agricultural changes (Jeffer- 
ies 2000; Fox et al. 2005; Van Eerden et al. 2005) on the 
wintering grounds. Concurrently, reproductive rates have 
dropped in some of these species, suggesting a density- 
dependent response in breeding output (Madsen, Cracknell 
& Fox 1999). Density dependence determines to what 
numbers these populations will ultimately rise, and is 
therefore also of management interest, because conflicts 
with agriculture on the wintering grounds have greatly 
intensified (van Eerden 1990). 

Voles and lemmings have long been known to exhibit 
multi-annual population cycles, especially in more north- 
erly regions (Elton 1942). Lemmings (Lemmus spp. and 
Dicrostonyx spp.) generally show population cycles with 
peaks every 3-5 years (Kokorev & Kuksov 2002; Gilg, 
Hanski & Sittler 2003; Krebs 2010). They play an impor- 
tant role in the Arctic food web, and are the preferred 
prey of generalist predators like Arctic foxes (Vulpes lag- 
opus L.), skuas and gulls (Lari) (Gauthier et al. 2004). 
Predation pressure by these generalist predators on alter- 
native prey like Arctic-nesting birds varies with lemming 
abundance, and the birds' reproductive output varies 
accordingly (Summers 1986). During the last 25 years or 
so, vole and lemming cycles are fading out in some north- 
erly regions, which has been attributed to changes in win- 
ter conditions (Bierman et al. 2006; Ims, Henden & 
Killengreen 2008; Kausrud et al. 2008; Gilg, Sittler & 
Hanski 2009). Through increased predation pressure on 
alternative prey, such collapses of lemming cycles have 
been suggested to have negative consequences for resident 
breeding birds (Kausrud et al. 2008). Most birds in areas 
where lemmings occur are, however, migratory (Gilg & 
Yoccoz 2010), and hence the question arises whether such 
negative consequences also exist for migratory birds. 

Dark-bellied brent geese (Branta bernicla bernicla L.; 
hereafter: brent geese) breed on coastal arctic tundra 
and winter in the temperate zone of western Europe 
(Fig. 1). This population has undergone a 20-fold 
increase since the 1950s (Fig. 2, top panel). Unlike many 




Fig. 1. Distribution of dark-bellied brent geese during winter, 
migration and breeding. Ovals indicate spring departure site in 
the Wadden Sea and the artic tundra in Western Taimyr, with 
the Pyasina river delta (with Mys Vostoclmy) located in the cen- 
tre (after Madsen, Cracknell & Fox 1999). 

other goose species, the population decreased after the 
mid-1990s, suggesting overcompensating density depen- 
dence (Ebbinge et al. 2002). However, its reproduction 
strongly fluctuates in concert with the approximately 
3-year lemming cycle (Blomqvist et al. 2002; Ebbinge & 
Spaans 2002; Summers 1986). When lemmings are 
scarce, reproductive success is invariably low, and there- 
fore a density-dependent relationship of the number of 
first-winter brent geese was only borne out after select- 
ing only those years with presumed low predation at the 
breeding grounds (Ebbinge et al. 2002). In contrast, an 
earlier analysis of the same population, only covering 
the growth phase of the population, concluded that 
there was no evidence for density dependence (Summers 
& Underhill 1991). Instead, these authors suggested that 
population growth could be predicted from a constant 
survival rate and a reproductive rate, oscillating with 
lemming abundance. 

However, even in years of high lemming abundance, 
breeding success of brent geese may be low (e.g. 1961, 
1967, 1994), so for a full understanding other factors like 
weather conditions should be considered as well (Ebbinge 
1989). A late spring in the arctic reduces Brent goose 
reproduction (Barry 1962; Lindberg, Sedinger & Flint 
1997; Spaans et al. 1998), and in summer, the gosling 
phase is a critical period in the survival of juveniles 
(Morrissette et al. 2010). 

To complicate matters further, pre-breeding conditions, 
and particularly those at the spring departure site may 
have carry-over effects on the breeding output of migra- 
tory birds (Ebbinge & Spaans 1995; Schmutz, Hobson & 
Morse 2006; Spaans et al. 2007). In the last few decades, 
these conditions have changed quite considerably, espe- 
cially on the main spring departure site in the Wadden 
Sea (Fig. 1; van der Graaf et al. 2006). Moreover, the 
Russian/Baltic barnacle geese (Branta leucopsis Bechstein), 
that make use of the same departure site in spring, have 
increased exponentially since 1960 (Fig. 2). These birds 
have also prolonged their spring staging period: whereas 
they used to leave at the end of March, many now stay 
until early May (Eichhorn et al. 2009). Competition with 
barnacle geese may reduce the conditions for spring 
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Fig. 2. Population development of dark-bellied brent goose and 
Russian/Baltic barnacle goose between 1955 and 2008. Lines indi- 
cate significant trends (see Results). Over this time period, brent 
goose breeding success as well as lemming abundance in Taimyr 
decreased, whereas spring temperature sums (expressed as grow- 
ing degree days), both at the spring staging site in the Wadden 
Sea (GDD n ) and at the Taimyr breeding grounds (GDD h ), 
increased. Average temperature during the gosling phase in sum- 
mer (T s ) did not change. The shaded area indicates the period 
since 1988 when overall lemming abundance and climate changed 
most markedly. 

fuelling of brent geese, because brent geese leave later, 
and are less inclined to graze on previously grazed patches 
than barnacle geese (Engelmoer et al. 2001; Stahl et al. 
2006). 

Hence, apart from population density, changes in sum- 
mer and winter season in the Arctic, as well as changes in 
winter and spring conditions outside the Arctic may all 
impact the performance of these arctic-nesting migrants. 



Here we apply a longer time series than Summers & Un- 
derbill (1991), including a phase of decline, to test 
whether their conclusion still holds that the population 
trajectory of dark-bellied brent geese can be understood 
from a reproductive rate that is solely dependent on lem- 
ming abundance, without invoking density dependence. 
Their hypothesis would predict that the population 
decline is due to faltering lemming cycles. 

Materials and methods 

STUDY POPULATION 

Dark-bellied brent geese breed mainly on the arctic tundra coast 
of Taimyr peninsula (Siberia), and winter 5000 km to the south- 
west along the western European coast. In spring, virtually the 
whole population gathers in the Wadden Sea for pre-migratory 
fuelling (Fig. 1). 

Population numbers were derived from the mid-winter (i.e. 
January) counts conducted under the auspices of the Goose 
Specialist Group of Wetlands International. In winter, this popu- 
lation confined to coastal areas in the Netherlands (20%), France 
(40%) and England (40%), with much smaller numbers in 
Denmark and Germany, and no major change in its winter distri- 
bution has occurred while the population size increased (Madsen, 
Cracknell & Fox 1999). Its concentrated occurrence enables full 
coverage of census counts, with accuracy likely to be within 10% 
of the true total (Fox et al. 2010). This is substantiated by the 
finding that annual survival estimated from these counts is within 
1% of that estimated from re-sightings of individually marked 
brent geese using a Jolly-Seber approach (Ebbinge 1992; Ebbinge 
et al. 2002). 

There is some evidence for local shifts in population distribu- 
tion in spring with changes in population size (Ebbinge 1992; 
Engelmoer et al. 2001). Such information from the breeding 
grounds is lacking, but given the loose colony structure and 
nomadic nature of brent geese in summer (Ebbinge 2004), we 
expect breeding density and hence intraspecific competition to 
increase with population size. 



FROM BREEDING SUCCESS TO POPULATION 
DEVELOPMENT 

We used an information-theoretical approach to determine the 
main environmental factors affecting population productivity 
(Morrissette et al. 2010). Then we used a general linear model 
containing the best explanatory factors to predict breeding suc- 
cess in a given year. We subsequently combined these with pub- 
lished estimates of survival to project population development 
under a few selected ecological scenarios in order to separate the 
effects of changes in lemming cycles and population density. 



SURVIVAL AND BREEDING SUCCESS 

Proportions of juveniles were assessed by ageing birds at all 
major wintering sites, with sampling effort (25-50% of the popu- 
lation) increasing with the increase in population. We used these 
proportions of juveniles, weighted for population size per coun- 
try, to subdivide the census population in different age classes, 
and to calculate population breeding success. Because brent geese 
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are capable of breeding at 2 years of age (Ebbinge et al. 2002), 
we subdivided the census N c in a given year / in three age classes, 
namely first-winter birds Ny, second-winter birds N s , and older 
birds N a . The number of first-winter birds was calculated from 
the proportion of first-winter birds p, as N/ it = p,. N Cit . The num- 
ber of second-winter birds was derived from the number of first- 
winter birds in the previous year as N S: , = s t . Nyt.i, where s is 
survival of birds > 0-5 year old (survival after the first 6 months 
was considered constant; Ebbinge et al. 2002). Survival was 
modelled to be lower in the first decades of the study as hunting 
was still allowed in Denmark until 1972 (Ebbinge 1985). Crude 
survival analyses yielded estimates of s of 0-81 before 1972, and 
0-85 after 1972 (Summers & Underbill 1991). Considering the 
whole period (both before and after 1972), s was, however, esti- 
mated to be higher (0-86) (Ebbinge et al. 2002; Summers & Un- 
derbill 1991), so we used s of 0-81 and 0-86 for before and after 
1972 respectively. Subtracting the first- and second- winter birds 
from the census number gave the number of older birds N a ,: 

N aJ = N CJ - (N f j + N s>l ) eqn 1 

Breeding success b, was subsequently calculated as: 

, N fj P< N C< n. 

b, = = ; r eqn 2 

b, was square-root transformed to normalize the data. 

LEMMING ABUNDANCE 

In brent geese, reproduction strongly fluctuates in concert with 
the lemming cycle (Blomqvist et al. 2002; Ebbinge & Spaans 
2002; Summers 1986), Relative lemming abundance (L) was 
derived from numbers of lemmings caught in snap-traps ordered 
along transects in different tundra habitats in summer (June- 
August). The main species is Siberian lemming (Lemmus sibiricus 
Kerr), with collared lemming (Dicrostonyx torquatus Pallas) being 
about an order of magnitude less abundant. L was classified into 
five categories according to the total number of lemmings caught 
per 100 trap-days (1: < 1; 2: 1-3; 3: 4-10; 4: 11-30; 5: > 30) 
(Kokorev & Kuksov 2002). We used the lemming index for arctic 
tundra in Western Taimyr (zone of c. 100 km along the coast of 
the Kara Sea between 80 and 90°E) (Chernov 1985) as provided for 
1960-2001 by Kokorev & Kuksov (2002). We complemented this 
series for 2002-2008 from snap-trap data collected at Mys Vost- 
ochny (74"08' N, 86°44' E) in the Pyasina delta, also in the arctic 
tundra of Western Taimyr (1200-1500 trap-days per summer) 
(Rykhlikova & Popov 2000). Due to its logarithmic nature, the 
lemming index is robust to slight changes in methodology: the lem- 
ming index determined at Mys Vostochny in 1993-1995 was identi- 
cal to that given by Kokorev & Kuksov (2002). No trapping data 
were available for 2003, but observations (http://www.arcticbirds. 
net) suggest a lemming index of 1, in accordance with the notion 
that the previous year, 2002, had been a lemming peak year (that 
did not grow to its full potential; Tulp 2007). 



CLIMATIC EFFECTS 

Pre-breeding conditions, and particularly those at the spring 
departure site may have carry-over effects on the breeding output 
of brent geese (Ebbinge & Spaans 1995; Spaans et al. 2007). In 



another Arctic-nesting goose species, body stores were correlated 
with the growing degree days (GDD) around spring departure 
(Duriez et al. 2009). As a proxy for departure conditions from 
the pre-migratory fuelling site we therefore took growing degree 
days (GDD lr ; subscript refers to wintering grounds). GDD W was 
calculated as the sum of average daily temperature above 0 °C 
measured at Leeuwarden weather station (53°12' N, 05°47' E) 
until 23 May, the mean date of mass departure of brent geese 
from the Wadden Sea. 

Brent goose reproduction is dependent on spring conditions 
(Barry 1962; Lindberg, Sedinger & Flint 1997; Spaans et al. 
1998), and because plant development is closely related to a 
spring temperature sum (Wang 1960), we used growing degree 
days in Taimyr (GDD/,; subscript refers to breeding grounds) as a 
proxy for arrival conditions in the breeding grounds. GDD h was 
calculated as the sum of average daily temperature above 0 °C 
measured at Dikson weather station (73°31' N, 80°20' E) until 19 
June, the mean date at which 50% of the nesting brent geese are 
present in Taimyr (Spaans et al. 2007). This measure was a good 
proxy for the date of snow melt: GDDh was negatively correlated 
with the date at which 50% of the snow at our transect in Mys 
Vostochny had melted (on average 20 June; Pearson's R = —0-75, 
N = 9, P u < 0 01). 

Mean temperature during the gosling phase is correlated with 
the annual productivity of another Arctic-nesting goose species 
(Morrissette et al. 2010). As a proxy for the conditions during the 
first month after hatching, we took the average daily temperature 
(T s ; subscript refers to summer) measured at Dikson weather sta- 
tion between 21 July (mean hatching date, Spaans et al. 2007) and 
20 August, when all surviving goslings have fledged. 

INTERSPECIFIC COMPETITION 

The Wadden Sea spring staging site is also used by barnacle geese 
during their preparation for migration to the Baltic Sea and the 
Barentsz Sea (Russia) (Eichhorn et al. 2009), and competition 
with barnacle geese may reduce the conditions for spring fuelling 
of brent geese (Engelmoer et al. 2001; Stahl et al. 2006). The 
numbers of barnacle geese M were also derived from mid-winter 
(i.e. January) counts (Ganter et al. 1999; Van der Jeugd et al. 
2009). Missing data in the early period with little population 
growth (10 years up to 1972) were linearly interpolated. 

DATA ANALYSIS 

For the period 1960-2008, we ran a general linear model with 
square-root transformed breeding success (Vft,) as the dependent 
variable, lemming abundance (L) as a categorical predictor, and 
as continuous predictors the natural log-transformed brent and 
barnacle goose population sizes of the preceding winter \i\(N,_i) 
and ln(M,_ ; ), respectively, the growing degree days in the 
Wadden Sea and in Taimyr (GDD„. and \og(GDD h ) respectively), 
and the average summer temperature (T s ). Because ln(7V,_/) and 
ln(M,_ ; ) were highly correlated (Pearson's R = 0-865, N = 49, 
P < 0 001), we did not consider these together to avoid multicol- 
linearity. Weather data were not complete for 1971 and 1972, so 
we used averages over 1960-2008 (\og(GDD b ) = 0-764 for 1972, 
and T s = 5-4 °C for 1971 and 1972). All possible 48 models were 
ranked according to Akaike's Information Criterion, corrected 
for small sample size (AICc; Burnham & Anderson 2002), using 
the glmulti-package in R (Calcagno & de Mazancourt 2010). 
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Akaike weights were calculated by summing over all models; 
models with AAICc > 2 were considered to have less than sub- 
stantial empirical support (Burnham & Anderson 2002). 

PROJECTED POPULATION DEVELOPMENT 

We used predictive modelling to separate the effects of ecological 
drivers, through their effects on breeding success, on population 
development. We used published estimates of annual survival s 
(see Survival and breeding success) and our modelled breeding 
success ft, to derive the expected number of birds N: 



Nfj = P,N aJ 
where 

As (see Survival and breeding success): 

N s ,i =s t Nf,t-i 

this gives an expected census number: 

Na = Nj„ + N s ,i + N a<l 



eqn 3 



eqn 4 



eqn 5 



eqn 6 



Such projections suffer from error propagation, and there- 
fore we did not model the whole time series. Because we 
were specifically interested in the factors causing the 
population to level off and even decline, we used the 
census with maximum numbers (1991) as the start for 
these projections. 

Results 

Since 1955, population size of the dark-bellied brent goose 
increased (ln(iV,), Pearson's R = 0-97, N = 55, P < 0-001), 
but clearly levelled off and even decreased in the 1990s 
(Fig. 2; year 2 was significant in polynomial regression: 
t S 2 = — 9-41 , P < 0-001). Concurrently with the increasing 
population size, brent goose' breeding success dropped 
(■4b „ Pearson's R = -0-29, N = 54, P < 0-05). The lem- 
ming abundance also showed a decrease, because peaks 
occurred less regularly, and were less pronounced (Fig. 2; 
Spearman p = -0-32, TV = 49, P < 0-05). In contrast to 



the brent goose, the Russian/Baltic barnacle goose popu- 
lation continued to grow in an exponential fashion 
(Fig. 2; (ln(M,), Pearson's R = 0-98, N = 40, P < 0-001). 
Spring temperature sums (i.e. growing degree days) at the 
breeding grounds in Taimyr increased (\og(GDD h ), Pear- 
son's R = 0-34, N = 48, P < 0-02), but average tempera- 
tures during the gosling phase did not change (T s ; 
Pearson's R = 0-17, N = 52, P = 0-23). Temperature sums 
at the spring departure site in the Wadden Sea also 
increased (GDD„.; Pearson's R = 0-538, N = 54, 
P < 0-001). The observed changes in lemming abundance 
and climate were all most prominent since 1988 (shaded 
area in Fig. 2). 

The most parsimonious model explained breeding suc- 
cess {-4b,) of brent geese with lemming abundance (L), 
brent goose population size, ln(JV,_/), and the growing 
degree days in Taimyr (GDD/,). This model clearly per- 
formed better (AAICc > 2) than all other models 
(Table 1). Considering all 48 models, these variables were 
also best explaining the observed variation in breeding 
success (Table 2). The size of the brent goose population 
itself featured much more prominently than that of its 
competitor, the barnacle goose, with a sum of Akaike's 
weights of 0-16 for ln(M,./) compared to 0-61 for ln(/V M ). 
(Note that w/s of ln(/V,_,) and ln(M M ) were slightly 
underestimated because these variables only featured in 
16 compared to 24 models for the other independent 
variables; ignoring multicollinearity and hence considering 
24 models for ln(N,_i) and ln(M,_i) too, their ir,- rose to 
0-66 and 0-27 respectively). Weights of evidence for effects 
of summer temperature at the breeding grounds and 
growing degree days at the spring staging site were also 
considerably lower (Table 2). Another way to test whether 
there was a density-dependent effect on breeding success 
is to perform a type three likelihood ratio test (Sokal & 
Rohlf 1995) for the best performing model. This revealed 
that the increment in log-likelihood by adding brent 
goose population size ln(/V,_/) to a model containing 
lemming abundance (L) and the growing degree days in 
Taimyr (GDD/,),vjas marginally significant (%i = 3-82, 
P ^ 0-05). So, breeding success of the brent geese was 
higher when lemmings were abundant, whereas it was 
lower in late springs (as indicated by a low GGD), and 



Table 1. Top 5 of explanatory models of breeding success (Vi>,) of dark-bellied brent geese according to Akaike's information criterion 
corrected for small sample size (AICc; n = 49). Each model is defined by the listed independent variables: lemming abundance L, spring 
conditions at the breeding grounds in Taimyr (expressed in growing degree days, GDD/,), spring conditions at the departure site in the 
Wadden Sea (GDD„), average summer temperature in Taimyr (T s ), and population size of competitor species (barnacle goose, ln(M,_ ; ), 
Models with a difference of AICc less than two are indicated in bold; K is the number of parameters, Ifflrn^x) is likelihood of model /' 
given the data, and w-, is its Akaike's weight. RMSD is the root mean squared deviation of observed on predicted values (of V/;,), and R 2 
the proportion of the variation in observed values explained by the predicted values (of 4b,) 



Model 


K 


Scaled deviance 


AAICc 


r /j(mi\x) 


»'/ 


Evidence ratio 


RMSD 


R 2 


(1) L - \n(N,_,) + \og(GDD„) 


8 


1 17 


000 


1 00 


0-35 


10 


0-22 


063 


(2) L - ln(W,.,) + \og(GDD h ) - T s 


9 


1-20 


2-33 


0-31 


011 


3-2 


0-21 


0-64 


(3) L - \n(N,.,) + \og(GDD h ) - GDD W 


9 


1-20 


2-61 


0-27 


010 


3-7 


0-22 


0-63 


(4) L + \og{GDD h ) - ln(Af w ) 


8 


117 


2-64 


0-27 


009 


3-8 


0-22 


0-62 


(5) L + \og(GDD h ) 


7 


114 


2-67 


0-26 


009 


3-8 


0-23 


0-60 
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Table 2. Weight of evidence for effect of the six test variables on 
reproductive output of dark-bellied brent geese (Vfi,). Standard- 
ized coefficients of the full model give an indication of their rela- 
tive effect size. Akaike weights are summed over all 48 models, 
and give an indication of how much support for an effect is pro- 
vided by the data. The Pearson's correlation coefficient indicates 
the simple correlation with yb t 



Standardized 



Variable 

Lemming class (L)* 
Growing degree days 

Taimyr, \o%(GDD h ) 
Brent goose population 

size, ln(N t _i) 

Summer temperature (T s ) 
Growing degree days 

Wadden Sea (GDD„) 
Barnacle goose population 

size, ln(M,.;) 



Pearson's R 


coefficients 


Wi 


0-72 


0-58 


100 


0-26 


0-29 


0-94 


-0-35 


-034 


0-61 


-0-24 


-013 


0-24 


-0-29 


—011 


0-24 


-0-31 


0-21 


016 



*these estimates were obtained by reclassifying lemming class as 
0 for L < 2 and 1 for L > 3. 

tended to be lower at large brent goose population sizes 
(Fig. 3). 

We used this most parsimonious model (model 1 in 
Table 1; for parameter estimations, see Table SI) to predict 
breeding success (from the environmental variables L, log 
(GDD h ) and ln(N t _])), and to project population develop- 
ment. The model nicely captured the decrease in the brent 
goose population after 1991 (Fig. 4, ml). Interestingly, the 
model without ln(A f ,. / ), i.e. excluding negative density 
dependence on breeding output (model 5 in Table 1), pre- 
dicted a levelling off, but not a decrease, in population size 
(Fig. 4). We investigated to what extent the decrease was 
due to the observed changes in lemming abundance L, 
which changed quite abruptly in 1988 (Fig. 2). Hence, we 
substituted the values for the 20 years including and follow- 
ing 1988 with those of the 20 years preceding 1988 in 
model 1. If lemming abundance had not changed after 
1988, brent goose population size would eventually have 
increased further (Fig. 4, mlc; note that this population 
increase was projected even though negative density depen- 
dence on breeding output was accounted for). 





Fig. 3. Breeding success (V/j) as a function of lemming abun- 
dance (L) and spring conditions expressed as growing degree days 
(°C) at the breeding grounds in Taimyr (logGDD h ) (a relatively 
low GDD represents a late spring). Breeding success increases 
with lemming abundance, and is higher when spring in Taimyr is 
early and, in general, when population size is low. (a) Observed 
values shaded by the brent goose population size \n{N t _j), increas- 
ing from white (< 100 000) to black (> 200 000). (b) Predicted 
values according to most parsimonious model (model 1 in 
Table 1) for population size of 15 000 (upper mesh) and 300 000 
(lower mesh). 



Discussion 

We have identified lemming abundance, the onset of 
spring at the breeding grounds and population density to 
be the best explaining factors of population productivity 
in dark-bellied brent geese over the last half-century. 
There was little evidence for carry-over effects arising 
from conditions at the main spring staging site. 

FALTERING LEMMING CYCLES 

The levelling off and even decrease in population size of 
brent geese is due to a series of years with low breeding 
output, which is mainly the result of faltering lemming 



cycles during the last two decades. Although these geese 
do not always produce many offspring when lemmings 
are abundant, the converse holds: in lemming trough 
years following a lemming peak year, when the generalist 
predators are numerous and are relying on alternative 
prey like brent geese eggs, breeding success is consistently 
low (Ebbinge & Spaans 2002; Gauthier et al. 2004). Our 
projective modelling suggests that if lemming abundance 
had not changed, brent geese population size would even- 
tually have increased further, despite density dependence 
operating. 

Changes in winter climate have been suggested to cause 
the lemming cycles to falter (Ims, Henden & Killengreen 
2008; Kausrud et al. 2008; Gilg, Sittler & Hanski 2009). 
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Fig. 4. Population development of dark-bellied brent goose 
according to the census (dots), and projected from 1991 (when the 
population size was maximal) onwards (black lines). Projections 
are based on crude survival estimates (see Methods) and breeding 
success predicted from lemming abundance and spring temperature 
sums at the breeding grounds (m5, model 5 in Table 1), or, 
additionally, the brent goose population size (ml, model 1 in 
Table 1). mlc = model 1 but with lemming cycles continuing as 
before 1988. For parameter estimates see supplementary Table SI . 

Lemmings can reproduce under snow cover, which pro- 
vides protection from predators and thermal insulation, 
while food plants are still accessible (Stenseth & Ims 1993; 
Millar 2001; Korslund & Steen 2006; Reid et al. 2012). 
Snow conditions are therefore a crucial determinant of 
lemming peaks (Ims, Yoccoz & Killengreen 2011). There 
are indications that these are changing in Taimyr. The 
snow cover period tends to be shorter, but snow depth 
has increased considerably in 1966-2007 (Bulygina, 
Razuvaev & Korshunova 2009), and hence these changes 
are not all negative for lemmings. In some recent years 
(2002 and 2008), lemming peaks did not seem to reach 
full potential, despite lemming nests being about as abun- 
dant under the snow (in 2008) as in a year that featured a 
lemming peak (2005) (Feige et al. 2012). It seems there- 
fore particularly relevant to investigate whether rain-on- 
snow or melt-refreeze events occur earlier or more often 
(Ye, Yang & Robinson 2008; Rennert et al. 2009), as 
these events may mean that more water is perforating 
through the snow, potentially jeopardizing the insulation 
of the lemmings, and perhaps ice layers on the tundra sur- 
face render feeding grounds inaccessible (Korslund & 
Steen 2006; de Raad, Mazurov & Ebbinge 2011). Such 
icing is also thought to be responsible for the lengthening 
of the lemming cycle from 5-8 years at Wrangel Island 
(Menyushina et al. 2012). 

COPING WITH LATE SPRINGS 

In the Arctic, the timing of snow melt, and hence the 
length of the summer season, varies considerably among 
years. Arctic-nesting geese start nesting at a somewhat 



later date in years when snow melts late (Prop & De Vries 
1993; Cooke, Rockwell & Lank 1995; Madsen et al. 
2006), but, relative to snow melt, commence nesting early 
in late springs (Barry 1962; Lindberg, Sedinger & Flint 
1997; Bety, Gauthier & Giroux 2003). Also, like in other 
bird species (Murphy 1986; Perrins & McCleery 1989), 
clutch size in geese is generally smaller in late springs 
(Barry 1962; Raveling 1978; but see Lindberg, Sedinger & 
Flint 1997; Madsen et al. 2006). Geese probably start 
nesting early relative to snow melt in late springs to 
enhance the survival of the goslings (Prop & De Vries 
1993). 

In their decision when to commence nesting, the birds 
face another trade-off, namely that between current and 
future reproductive success (Daan, Dijkstra & Tinbergen 
1990). In many geese and duck species, only females 
incubate the eggs, and the body weight at the end of the 
incubation is supposedly near the critical boundary (Drent 
et al. 2003). Laying a smaller clutch may compensate for 
the lower pre-laying condition of the female, leading to 
equal body weights at the start of incubation irrespective 
of the onset of spring (Se'ne'chal, Bety & Gilchrist 2011; 
Spaans et al. 2007). Such compensation may be only 
partial, as suggested by the finding that female black 
brants {Branta b. nigricans Lawrence) weighed less 
during incubation in a later spring (Eichholz & Sedinger 
1999). 

Commencing nesting at a higher snow cover in late 
springs would mean that poorer feeding opportunities 
may extend well into the incubation phase (Prop & De 
Vries 1993). Then female geese would have to compensate 
for poorer feeding conditions in the early incubation 
phase by spending more time feeding (black brant: 
Eichholz & Sedinger 1999; dark-bellied brent geese: 
B.A. Nolet, R.A. Bom, J. de Fouw, P.P. de Vries, B.S. 
Ebbinge, unpublished). Limiting the loss in body condi- 
tion would clearly enhance the female's survival prospects 
and future reproductive success. Because leaving the nest 
by the female is risky for dark-bellied brent geese because 
of egg predation (Spaans et al. 2007), egg predation is 
expected to be heavy in late spring. So, egg survival seems 
to be traded-off against both gosling and female survival 
to make the best of a bad job. 

All in all, reproduction is severely reduced in late 
springs. Such late springs are predicted to become less fre- 
quent (Tulp & Schekkerman 2008), and this may partly 
compensate the effects of faltering lemming cycles (see 
above). 

DENSITY DEPENDENCE 

111 an earlier study, Summers & Underhill (1991) did not 
find evidence for density dependence operating during the 
growth phase of the dark-bellied brent goose population 
(1955-1988). At the end of that period population size 
was just above 200 000 individuals. Our analysis, which 
covers another 20 years in which the population roughly 
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fluctuated between 200 000 and 300 000 individuals, 
shows that there was a negative, albeit weak, effect of 
population size on breeding output. This suggests that 
density dependence only starts to have an effect at a pop- 
ulation size above 200 000 individuals. 

Recent analyses suggest density dependence is acting in 
some Arctic-nesting goose species like the Svalbard barna- 
cle goose (Trinder, Hassell & Votier 2009), but not in oth- 
ers like the greater snow goose (Anser caerulescens 
atlanticus Kennard) (Morrissette et al. 2010). During the 
growth phase of the dark-bellied brent goose population, 
an increasing proportion of the birds staged at a less pre- 
ferred site, suggesting that the preferred site was filled to 
capacity (Ebbinge 1992). If, as a consequence, a larger 
proportion of females leaves with a relatively low condi- 
tion, this might affect subsequent population productivity 
(Ebbinge & Spaans 1995). However, in our analysis we 
did not find evidence for such carry-over effects (see 
Carry-over effects), and we think that density dependence, 
noticeable on the level of the entire population, acts 
mainly on the summer grounds. 

The main nesting habitat of brent geese are coastal 
islands, that are usually also frequented by Taimyr gulls 
(Lams fuscus taimyrensis Buturlin). Only in lemming peak 
years, considerable numbers also nest on the mainland 
around nests of snowy owls (Bubo scandiacus L.) that 
only breed in such years (Underhill et al. 1993; Summers 
et al. 1994). There are indications that breeding distribu- 
tion is despotic in both nesting habitats, with clutch sizes 
being larger close to gulls' nests (at least in predation 
years; J. de Fouw et ah, unpublished) and close to snowy 
owls' nests (van Kleef et al. 2007). The negative density 
dependence may arise from these best nesting territories 
being filled to capacity (Rodenhouse, Sherry & Holmes 
1997). 

CARRY-OVER AND INTERACTION EFFECTS 

In a classic study, it was shown that dark-bellied brent 
geese that are heavy upon departure from the Wadden 
Sea spring staging area have a higher probability to breed 
successfully (Ebbinge & Spaans 1995). Contrary to the 
notion of 'Arctic amplification' (Solomon et al. 2007), the 
climatic changes on this spring staging site in the temper- 
ate region were even more pronounced than those at the 
Arctic breeding grounds. However, measured across the 
whole population we did not find an effect of spring 
weather conditions in the Wadden Sea on reproductive 
output. We also did not find evidence for a depressing 
effect of the increase in the competitor species, the barna- 
cle goose. Competition between barnacle and brent geese 
is intense on the departure site in the Wadden Sea, with 
brent geese suffering more from barnacle geese than the 
other way around (Stahl et al. 2006). However, these spe- 
cies are spatially segregated at the breeding grounds 
(Cramp & Simmons 1977). Brent goose and barnacle 
goose population development were partly running 



parallel, and hence the effect of barnacle geese may have 
been masked by even stronger effects of brent geese popu- 
lation density itself. 

To avoid possible spurious effects we had to reduce the 
number of models under consideration, and we did not 
investigate possible interaction effects. For instance, the 
effect of a late spring at the breeding grounds might be 
less detrimental if it coincides with a late spring at the 
departure site. In general, a late spring at the departure 
site may lead to a later departure (Duriez et al. 2009; 
Stirnemann et al. 2012), reducing the costs of waiting at 
the breeding grounds. Such effects may be small in 
dark-bellied brent geese, however, because departure is 
probably more triggered by wind conditions (Green 
2001). 

Conclusion 

Lemmings have long been known to affect, indirectly, 
population productivity of Arctic-nesting migratory birds 
(Summers 1986). More recently, population dynamics of 
resident bird species were shown to change in concert 
with changes in lemming dynamics, but this might be due 
to direct effects of snow hardness on both lemmings and 
resident grouse species (Kausrud et al. 2008). Our analysis 
points to faltering lemming cycles, perhaps caused by 
changes in the Arctic winter, as the main factor causing 
changes in both population productivity and population 
size of a migratory bird. Although in the last decade our 
knowledge about the dynamics of these rodents has 
improved, different factors seem to play a role at different 
places (Gilg, Hanski & Sittler 2003; Kausrud et al. 2008; 
Krebs 2010; Menyushina et al. 2012). Hence, much more 
research is needed before longer term predictions on lem- 
ming abundance and all its dependencies, sometimes 
apparently stretching 5000 km, can be made. 
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